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The evolution of the non-adiabatic pressure perturbation during inflation driven by two scalar 
fields is studied numerically for three different types of models. In the first model, the fields have 
standard kinetic terms. The other two models considered feature non-canonical kinetic terms; the 
first containing two fields which are coupled via their kinetic terms, and the second where one field 
has the standard kinetic term with the other field being a DBI field. We find that the evolution 
and the final amplitude of the non-adiabatic pressure perturbation depends strongly on the kinetic 
terms. 



I. INTRODUCTION 



Inflation is the most successful theory when describing primordial perturbations in the universe. These primordial 
density perturbations are generated from quantum fluctuations of the field which has driven inflation (or, as an 
extension, of the perturbations in the curvaton field). This theory is consistent with observations. Details about 
the nature of the primordial perturbations depend on the details of the theory, such as whether there was only one 
field responsible for inflation, the form of the potential and interactions between fields. Observations of the cosmic 
microwave background radiation (CMB) constrain the amount of isocurvature perturbations to high accuracy (of 
order of 10% of the curvature perturbations) and thereby ruling out a large number of inflationary models. 

It has been recently pointed out that one important quantity to be studied in detail is the non-adiabatic pressure 
perturbation. Apart from being an essential ingredient to determine the evolution of the curvature perturbation, 
it was shown to source vorticity perturbations at second order in perturbation theory [T]. To be more precise, the 
evolution of the second order vorticity uj2ij was shown to obey 
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where on the right hand side the first order quantities 8p\ (the energy density perturbation at first order) and 8P na d 
(the non-adiabatic pressure perturbation) appear. If the non-adiabatic pressure perturbation vanishes at first order, 
the vorticity decays not only at first order, but also at second order. However, if SP na d is non-zero, which is usually 
the case in multifield inflation, vorticity at second order is sourced by the non-adiabatic pressure perturbation. This 
opens up more possibilities of tests of inflation, as the non-zero vorticity can source -B-mode polarisation of the CMB 
photons. It is therefore important to understand the behaviour of <5P na d in diverse models. In addition to sourcing 
vorticity, the non-adiabatic perturbation affects also the evolution of the curvature perturbation on super-horizon 
scales, and thereby influencing the predictions for non-Gaussianity in these models, see e.g. [H [3] and references 
therein. 

In [4] , the evolution of the non-adiabatic pressure perturbation has been studied in detail for theories with canonical 
kinetic terms with different choices of potentials. In this paper, we consider the evolution of <5P na( j in models with 
non-canonical kinetic terms. In particular we consider a theory in which the kinetic terms are coupled as in [5HZ], 
and a model in which one field has the standard form with the other being a DBI field (see e.g. [51-HTj] and references 
therein) . 

The paper is organised as follows: In the next Section ( |TT| we write down the theories considered in this paper and 
derive the evolution equations for the perturbations. In Section|ITT]we present the results of our numerical calculations. 



Our conclusions are presented in Section IV 
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II. THE MODELS 



We will consider three models in this paper. The actions are given as follows: 
1. Two of the models we consider are described by actions of the form 

r [ Ji/f 2 i i 

S = I d 4 x^ -^R- ^d^d^- -e^g^d^d.X-V^x) 
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where in the first model &(</>) = and in the second model 

K4) = P4> 

and /3 is a constant. 

2. The third model we consider contains a scalar field with a canonical kinetic term and one DBI field. The action 
is given by 



S = J d A x^—g 
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7 is the warp factor describing the shape of the extra dimensions. Both A and fx are constants. 
In all the actions given, Mpi = (87rG) -1 / 2 is the reduced Planck mass and R is the Ricci scalar. 
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A. Background Equations of Motion 



We assume a spatially flat Friedmann-Robertson- Walker, FRW, spacetime 



ds 2 = -dt 2 + a 2 {t)5 i jdx i dx j 



where a(t) is the scale factor. 



(4) 



1. Kinetic Coupling 

The equations of motion for both fields and the Friedmann equations are given by [SJ [TT] 
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2b T 



X + (3H + 2b^) X + e~ M V x = 



and 
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where V<j, = dV/dcf> and similarly for V x , = db(cf>)/d(f> and the Hubble parameter is given as H = a/a. These 
equations reduce to the two standard field case when the kinetic coupling is set to zero. 
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2. DBI Field 

The equations of motion for the helds for this model are 

4> + 3H + V 4> = {) 
X + '3H~f- 2 x + ^ (1 - 3 7 ~ 2 + 2 7 - 3 ) + j- 3 V x = 

and Einstein's equations give 



3H 2 = ^ 2 + j( 1 -l) + V 
-2H = 2 + 7X 2 
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where f x = df jd\- The dots represent differentiation with respect to cosmic time t. 

B. First Order Perturbation Equations 

We turn now our attention to the first order perturbation equations and work in the longitudinal gauge, in which 
the line element is given by 



ds 2 = -(1 - 2^)dt 2 + a 2 (l + 2^)6 ij dx i dx j . 
First, we decompose the fields into their homogeneous and perturbed parts 

cf)(t, x) = 0(t) + 5<P(t, x), X (t, x) = x{t) + S X (t, x) 



(13) 
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and we will be working with the Fourier components of the perturbations, 8<pk(t) and 5xk(t) and will be omitting 
the subscript k to shorten the expressions. 



The perturbed Einstein field equations for the model concerning the scalar field and the DBI field are 
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We find the Klein-Gordon equations for the field perturbations 
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It is convenient to work with the gauge-invariant Mukhanov-Sasaki variables [HE3], defined by 
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The gauge-invariant form of the perturbation equations are 
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with the coefficients B$, B x , Cd,*, Ca, x , C xx , C x s, in the equations are as follows 
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We will later calculate the power spectrum of the curvature perturbation. The curvature perturbation is defined by 

* . (29) 
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The pressure perturbation SP is composed of adiabatic and non-adiabatic parts 

SP = c 2 s 5p + SP nad , 

P ; c 



(30) 



where c 2 = K is the adiabatic sound speed, Sp and 5P na d are the perturbations in the energy density and non-adiabatic 
pressure, respectively. 

For each model the adiabatic sound speed is given as follows: 



1. Two standard scalar fields 



c 2 = 1+ 2(¥ + y 
3H{<j> 2 + X 2 ) 



(31) 



2. One scalar field and one with a kinetic coupling 



c 2 = l 
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3. One scalar field and one DBI field 
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For the last two models, the sound speed reduces to the two standard scalar field model when j3 — and 7 = 1. 
Finally, the gauge-invariant entropy perturbation p~4l [15] is defined as 



(34) 



S = -pSP n ad- 



Due to the complexity of these equations, for all models considered, we will evaluate them numerically following 
the method outlined in [TO UMTS] , 
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III. RESULTS 



We will now describe the results of our numerical calculations. To be concrete, we consider all models presented in 
Section [n] with the double quadratic potential |19j which will be rewritten into the form 

V{<j ) ,x) = \ml{<j> 2 +T 2 x 2 ) (35) 

where 

r = — £ . 36 

As usual, to^ and m x are the masses of the two scalar fields. All the following plots are shown in the WMAP pivot 
scale where k = 0.002Mpc _1 [20] . 



1. Two Standard Scalar Ftelds 



In this model, we consider the case where T — 7.0 as studied in [TTJ[5T]. Further conditions were applied, = 
1.395 x 10~ 6 Mpi in order to match the numerical simulation to WMAP measurements of the curvature perturbation. 
Initial conditions for the two fields are xo — 0o = 12-0 as in [4]. 




FIG. 1: The background dynamics of the two fields, 0-field (purple) and x-field (blue), for the two standard scalar field case. 
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FIG. 2: A comparison of the curvature and entropy power spectra, Viz (blue, dashed) and Vs (purple) on the left, and the 
power spectra of the total and non-adiabatic pressure perturbations, Vsp (red, dashed) and Vsp nad (green) on the right. These 
are for the two standard scalar fields case. The dashed lines indicate horizon crossing. 



We see in Figure [T] that the x-field reaches the minimum of the potential, resulting in the (/>-field dominating for 
the last 30 e- folds of inflation. This is reflected in Figure [2] which at this point there is a rise in the power spectrum 



G 



of the curvature perturbation. The amplitude of the entropy perturbation S reduces dramatically when the exchange 
in field dominance occurs at 40 e-folds. It starts to gradually increase until a peak is reached at N = 70, at which a 
drop is experienced. At the end of inflation, we see the magnitude of the entropy perturbation is many orders smaller 
than the curvature perturbation, specifically S ~ 10~ 31 whereas 7Z ~ 10~ 9 . 

The behaviour of the non-adiabatic pressure and entropy perturbation power spectra generally agree at horizon 
crossing at N = 15 and onwards. 

We find that Figure [2] is in agreement with Figures 1 and 2 from [J]. 



2. One Standard Scalar Field and One Field containing a Kinetic Coupling 

We have two possible cases that will arise for the background; one in which the 0-field reaches the minimum of the 
potential well before the x-field, and vice versa. We will begin by considering the latter. 

In this scenario, the first 50 e-folds are dominated by the 0-field until at which point, there is an exchange in the 
field contributions, leaving the remaining ^-field until the end of inflation. 

The ratio of the field masses is T = 0.3 with = 6.395 x 10~ 6 . Initial conditions are \o — 12-0 and <fio — 11.0. In 
this case, the kinetic coupling is = 0.1. 
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FIG. 3: The background plot of the <^-field (purple) and x-field (blue) for the case where 0-field is initially dominant during 
the inflationary period. 
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FIG. 4: A comparison of various power spectra. Vn (blue, dashed) and Vs (purple) are displayed on the left and the right, 
Vsp (red, dashed) and Vsp nad (green). 



The behaviour of the entropy and non-adiabatic component of the pressure perturbations is significantly different 

However, 



to the model previously considered in Section pi 1 
before the end of inflation 



there is a rise and fall in S during the last 6 e-folds 
This feature can also be seen in Figure [2] Like the two scalar field model previously 
studied, the amplitude of S is many times smaller than 1Z. At the end of inflation, the final amplitudes for S and 
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Vsp nad ar e S ~ 10~ 26 and Vsp nad ~ 10~ 50 . We find that the the final entropy amplitude for this model is 10 5 larger 
than found in the model only considering canonical scalar fields. 

We now consider the other possible case where the x-field at first dominates the inflationary period. For this, the 
parameters used are V — 6.0 with = 1.005 x 10~ 6 . Starting conditions for the two fields are xo = 7.4 and 4>o = 7.5. 
With these parameters, the final 8 e-folds are dominated by the canonical scalar field. 
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FIG. 5: The background dynamics of the <p-&eld (purple) and x-held (blue). 
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FIG. 6: A comparison of the various power spectra, Vv, (blue, dashed) and Vs (purple) on the left, and Vsp (red, dashed) and 
Vsp nad (green) on the right. 



As expected, this choice of parameters has resulted in an one scalar field scenario at the end of inflation, similar 
to Section III 1 Due to this, the shape of the power spectra for all considered quantities will be similar. However, 
there is a slight difference in the last few e-folds in Vs (and Vsp nad ) which relates to the behaviour of the remaining 
canonical scalar field. S is seen rising and falling during the remaining e-folds of inflation in Figure [2] and Figure [4] 
Instead we see in Figure [6] that S continues to decrease until slow-roll is no longer satisfied. The final amplitudes for 
the entropy and non-adiabatic pressure perturbations are significantly larger than for the two standard scalar fields 
case. The values for the entropy and non-adiabatic pressure perturbation amplitudes are the same as those found in 
the previous case (where the x-field reaches the potential's minimum before the 0-field). 

There appears to be no difference in the curvature and entropy amplitudes between the two cases considered here 
in Section UTOI 



3. A Scalar Field and DBI Field 



As with the previous model containing the kinetic coupling, this DBI model will also be associated with the same 
two scenarios, one in which the DBI field decays before the scalar field and the reverse. For all the cases considered, 
the parameters that belong to the DBI model hold the following values, A = 2.0 x 10 12 and fi — 0.2 |S]. 
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First we consider the case where inflation is initially dominated by the scalar field 4>- In this scenario, T — 2.0 with 
TO^ = 115.6 x 10 -7 and the initial values for the two fields are \o = 2.0 and 4>q = 12.0. These parameters were chosen 
so the curvature perturbation amplitude agrees with WMAP measurements. 




FIG. 7: The background dynamics of both the 0-field (purple) and x-field (blue) are shown on the left. In this scenario, the 
DBI field is the first to fall to the minimum of the potential well. On the right shows the evolution of the parameter 7 during 
inflation. 




FIG. 8: A comparison of the power spectra; on the left are Viz (blue, dashed) and Vs (purple) and the right, Vgp (red, dashed) 
and Vsp nad (green). 



From Figure [7j there is an exchange in the field contributions at N — 44 and this is displayed in Figure [8] through 
the fall in the amplitude of the entropy perturbation. Unlike all the previous actions that were examined, the entropy 
and non-adiabatic pressure perturbation amplitudes for this particular model of DBI inflation, does not decrease 
during the final few e-folds of inflation. Furthermore, the final value of the entropy amplitude, S ~ 10~ 16 is markedly 
greater than found in Section |III 1| Similarly, this increase in amplitude is also found in the non-adiabatic pressure 
perturbation Vsp„ ad ~ 10~ 40 , whereas for the two scalar field case Vsp nad ~ 10 -54 . 

In the other case where the DBI field decays before the scalar field, parameter values are T = 35.1 where the mass 
of the </>-field is 65 x 10~ 8 . The fields have the same starting values as in the previous DBI field case. 
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FIG. 9: The background dynamics of both the 0-field (purple) and x~ ne ld (blue) are displayed on the left. In this case, the 
DBI field is the first to fall to the minimum of the potential well. The evolution of 7 is shown on the right. 



At first, the DBI field dominates the inflationary period until it reaches the minimum of the potential well and 



oscillates, at which the </>-field will become dominant. This is shown in Figure 10 through the drop in the amplitude 
of the entropy perturbation 12 e-folds before the end of inflation. 
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FIG. 10: A comparison of the power spectra. Viz (blue, dashed) and Vs (purple) are shown on the left and Vsp (red, dashed) 
and Vsp nad (green) are on the right. 



The behaviour of all power spectra are similar to those seen in the two canonical fields and kinetic coupling models, 
where the x-field is the first to reach the potential's minimum. This is due to the remaining 0-field experiencing 
slow-roll. As expected, the drop in the entropy power spectrum at N = 47 is due to the DBI field reaching the 
potential's minimum. Final amplitudes for the entropy and pressure (non-adiabatic) perturbations are Vs ~ 10 -37 
and V 5 p nai ~ 1(T 62 . 

Both these values differ greatly than those found for the previous case (when the 0-field is the first to fall to the 
minimum of the potential). Furthermore, the final amplitudes for S and SP na d for this scenario are much smaller 
than in the two scalar field case in Section fill 11 



IV. CONCLUSION 



In this paper we studied the evolution of the non-adiabatic pressure perturbation 5P na( j in multi-field inflation. We 
considered three models containing two massive scalar fields, but each model differs by the form of the kinetic terms. 
In all models, inflation is driven initially by the two fields. However, one of the field approaches the minimum of its 
potential before the end of inflation, so that inflation subsequently is driven purely by the second field. Our numerical 
results for the model with standard kinetic terms agree with [3]. Our main results can be summarised as follows: 
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• The nature of kinetic terms affect the evolution and the final value of <5P na d- In the case of the fields coupled by 
kinetic terms (model 2 with j3 non-zero) the amplitude of <5P na d does not depend on which of the field approaches 
zero first. However, the amplitude is roughly five orders of magnitude larger for (3 = 0.1 when compared to the 
canonical case (J3 = 0). 

• The evolution of Vs is affected by how fast the first field approaches the minimum of its potential. We generically 
observe an increase in Vs, which influences the results of V-jz- We find this effect to be largest in the theory 
with standard-kinetic terms. In all other models, the first field approaches zero much more gradually, resulting 
in only a slight increase in Vs ■ 

• In the third model, if the DBI field drives the last e-folds of inflation, we find that <5P na d is substantially larger 
than in other cases. On the other hand, if the DBI field is not significant in the later stages of inflation the 
results are comparable to the standard case, as expected. In this model, the amplitude of Vs is largest. 

The relevance of our results come from the fact that <5P n ad sources vorticity at second order, which affects predictions 
for the B-mode polarisation of the CMB. The model with the DBI field driving the last e-folds of inflation predicts the 
largest amplitude for <5-P na d at the end of inflation, which implies that this model predicts a larger source for vorticity 
than the other models. A comprehensive analysis of pre- and reheating in these models is necessary to estimate the 
amount of entropy perturbations in the radiation dominated epoch. The results will be model-dependent, because 
the details depend on the coupling of the inflaton-field(s) to the matter fields. 
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Appendix A: One Standard Scalar Field and One Field containing a Kinetic Coupling 

1. Perturbation Equations 



The perturbed Einstein field equations for this model are 



^ + + {2H + 3H 2 )^ 
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All these equations will reduce to those for the two canonical scalar fields when the kinetic coupling b((f>) is set to 
zero, in this case (3 = 0. 
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